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At nonzero temperature, pions propagate through a thermal medium at less than 
the speed of light. About low temperature, this effect begins not at ~ T 2 , but at 
next to leading order, ~ T 4 . We also derive the generalization of the relation of 
Gell-Mann, Oakes, and Renner to nonzero temperature. 



1 General analysis 

In this proceeding we give a pedagogical explanation of recent work of ours.0 
In this section we show how the velocity of pions, or more generally Goldstone 
bosons, changes with temperature in a thermal distribution. 

This in itself is pretty trivial, and familiar from other contexts. For exam- 
ples, the appropriate analogy to Goldstone bosons in a theory which is relativist 
tically invariant at zero temperature are spin waves in an antiferromagnet. cr@ 
This is because the dispersion relation for such spin waves is U) = vp, where 
lu is the frequency, p the magnitude of the spatial momentum, and so v is the 
velocity. (Spin waves in a ferromagnet behave as u> ~ p 2 .) In a relativistically 
invariant system, of course v = c, but in antiferromagnets, the velocity v is just 
some parameter that depends upon details of the couplings, etc. Thus for anti- 
ferromagnets, it is completely unremarkable that v changes with temperature, 
since the couplings change with temperature as well. 

This is in agreement with Goldstone's theorem, which tells us that for 
the Goldstone modes, the inverse pion propagator vanishes at zero spatial 
momentum, p — 0. □ Thus it is fine if ui vanishes like some constant v times p, 
where v changes with temperature. _ 

Less trivially, one can also understand the damping of pions or spin waves. □ 
Even without calculation, we expect that the damping rate vanishes at zero 
spatial momentum, since by Goldstone's theorem the inverse pion propaga- 
tor must vanish at zero momentum, for both the real and imaginary parts. 
How precisely the imaginary parts, and so the damping rate, vanishes is an 
interesting question which we face when we get into the gory details below. 

A similar problem is the nropagation of light in a medium with an index 
of refraction greater than one.Q The analogy is imprecise, though, because in a 
thermal bath Debye screening generates a finite correlation length for photons, 
whereas in the broken symmetric phase, Goldstone bosons always have an 
infinite correlation length. 
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A better example is the propagation of light in a background magnetic 
field, where at two loop order the velocity of light is less than thai in vacuum 
by an amount ~ a 2 B 2 /m 4 , where m is the mass of the electron. U There is a 
technical similarity here as well. The deviations from the speed of light are 
proportional to the energy density ~ B 2 for photons in a background B field, 
and to the free energy density ~ T 4 for thermal pions. For the former, this 
is unsurprising, since one wouldn't expect anything to depend upon B, only 
on B 2 . But for thermal pions, that there aren't terms pi ~ T 2 is at first sight 
unexpected. However old arguments in a different guiseQ tell us that to leading 
order, ~ T 2 , lorentz covariance is manifest in a pion gas, so v = 1. Indeed, it is 
amusing that for thermal pions, the deviation of v 2 from unity is proportional 
to 8tt 2 /45T 4 , Eqs. (HJ) and (||), which is 16/3 times the free energy density 
of a relativistic pion gas. 

Quite remarkably, this is analogous to the behavior of spin waves in ferro- 
magnets at low temperature, where the dispersion relation is uj(k) — c(T)k 2 . 
As first shown by DysontH, c(T) changes with temperature but only at next to 
leading order in an expansion about T = and the correction is proportional 
to the free energy density of the gas of spin waves. Presumably, the same 
is true for other systems with Goldstone modes at low temperature. Also, 
perhaps this a hint of the existence of a more general relation, which is valid 
nonperturbatively; analogous, say, to the expression for the speed of sound in 
a medium in terms of the pressure and density. 

Consequently, the basic physics which we consider is hardly revolutionary. 
Nevertheless, we think that we have a novel way of looking at these phenomena, 
for instead of considering inverse (pion) propagators, we use PC AC, and speak 
of pion decay constants. 

At zero temperature, the pion decay constant f n ~ 93 MeV is defined by 

<0K|7T b (P)) =ius ab p» , (1) 

where A% is the axial vector current, and the pion has momentum P^ — (p°,p). 

Once stated, our basic point is obvious. At nonzero temperature, because 
of the presence of the medium, there are, in general, two distinct pion decay 
constants. The timelike component of the current has one, 

(0\A°\* b (P)) T = ifi5 ab p° , 

and the spatial part of the current has another, 

(0\Ai\ir b (P)) T = if s J ah f • (2) 

At this point we should be careful and qualify exactly in which regime 
we are computing. First, all matrix elements are computed at a tempera- 
ture T in the imaginary time formalism; then, the timelike component of the 
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momentum, p° , is analytically continued from euclidean to Minkowski values, 
p° = —iu) + + . The two pion decay constants, fl and are denned about 
zero momentum, uj and p — > 0; typically, each has both a real and an imaginary 
part. 

In a colloquial fashion, we only deal with soft, "cool" pions. Soft means 
that all components of the momenta are small relative to (the real parts) of 
fl and "Cool" means that the temperature is well below that for the 
restoration of chiral symmetry. This is certainly a well defined regime, where 
the sigma meson is heavy and its abundance is Boltzmann suppressed. Once 
the sigma becomes light, with a mass less than the f^'s, then one enters the 
regime of the critical point, and our analysis does not apply. 

We certainly aren't the first people to realize that there are two pion 
decay constants at nonzero temperature: that's an old story in nonrelativistic 
systems; B for pions, it_was pointed out before by Kirchbach and Riska,l3 by 
Thorsson and Wirzba, B and probably by others as well. 

What we want to emphasize is that whenever there are two, distinct pion 
decay constants, then the propagation of pions must move off of the light cone. 
The argument is elementary: whether or not the vacuum respects the chiral 
symmetry, the symmetry is still good on the pion mass shell. Thus if we take 
the derivative of the axial current, we can use it to define the pion mass shell: 

ft Pa + /w P 2 — 0|tt mass shell ■ (3) 

This expression has both a real and an imaginary part. The mass shell condi- 
tion for pions is then, 

p° = -iu - 7 , (4) 

where the real part gives 

2 2 2 ^ ^ e fn 2 rr\ 

w=vp "rc7^- (5) 

and the imaginary part fixes 

1 



7 



^^(+W^ 2 -I-/> 2 )>0, (6) 



In doing so, we assume that the imaginary parts are small relative to the real 
parts, 

Im/^«Re/^. (7) 

One could probably gain some general understanding of what happens when 
this condition doesn't hold, and pions are strongly damped, but we haven't 
tried. 
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With Eq. (§) in hand, we can understand the possible behavior of the 
damping rate about zero momentum. The simplest possibility is that the 
imaginary parts of /* and are nonzero at zero momentum; then the damping 
rate vanishes linearly as p — > 0. One can show that in this case, the imaginary 
part of the pion self energy vanishes quadratically about zero momentum, so 
Goldstone's theorem is evidently obeyed. However, in some cases in which 
calculations have been done, the damping rate vanishes faster than linearly. 
For example, in the linear sigma model discussed in the next section, it vanishes 
exponentially. This is special to the linear sigma model at the order treated, 
since the only intermediate state through which the pion can scatter involves 
a virtual sigma meson, which is heavy, and so Boltzmann suppressed. Also, 
for spin waves in an antiferromagnet, the damping rate vanishes quadratically 
about zero momentum, cl _ 

We quote, without derivation, our results EJ on the generalization of the 
relation of Gell-Mann, Oakes, and Renner to nonzero temperature. Using 
methods of Shore and Veneziano,EJ we find that the dynamic pion mass, which 
is the pole for complex p° at p = 0, is 



a 2m(qq) T 
(Re/*)*" 



This is the same expression as at zero temperature, except that instead of f v 
entering, it is the real part of /* which appears. Of course, both Re f* and the 
value of the quark chiral condensate, (qq)T, are functions of temperature. A 
relation like this was first obtained by Thorsson and Wirzba, □ although they 
just wrote /* instead of its real part. 

It is worth emphasizing that this is the dynamic pion mass, and not the 
static pion mass. The latter is given by the pole for zero frequency, p° = 0, in 
the complex p plane. This mass is just 



m staUc = m* ^ (Q) 



V 



and so is less than when v < 1. This is interesting, because what is 
typically measured in numerical simulations of lattice gauge theory are the 
static and not the dynamic masses; the static mass is then an upper bound on 
the dynamic mass. 

2 Lowest order 

We now consider where these effects fixst appear io-an expansion about zero 
temperature. Using either a nonlinearLLj or a linearis sigma model, to leading 
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order in T 2 // 2 , fl(T) — f%(T), so pions move undamped at the speed of light. 
The reason for this was given by Dey, Eletsky, and Ioffe, Q who considered 
the thermal average of the two point function of either vector or axial vector 
currents. For completeness, we redo their analysis in terms of the thermal 
average of the one point function^], 

(0Kk> T - E>|e _* /T|n) do) 

At low temperature, T <C f n , the states \n) contain only pions. Further, to 
lowest order in T 2 // 2 we can truncate states with the fewest number of pions 
to obtain 



/d 3 k 1 
___ n(fc)(7rC(fc) |^| 7 r c (fc);^) T =o + ... 

where k = |k| and n(k) is the Bose-Einstein statistical distribution function. 
The Bose-Einstein distribution function is built up by summing over states: in 
the bra state, first one pion with momentum k, weighted as e -fc / T , then two 
pions with the same momentum, weighted as e~ 2k / T , etc. In bra states with 
more than one pion, all pions beyond the first are completely disconnected from 
the matrix element, and only contribute through their statistical weighting. 
Summing over all such pions then generates e~ fc / T /(l — e~ fc / T ) = n(k). Using 
Aq as an interpolating field for 7r a , and the canonical commutation relations 
of current algebra, the thermal pions n c (k) are eliminated to obtain 

$> c (fc)K|7r c (fc);7r b > T=0 ~ -^<0KK b > T=o . (12) 
To leading order in T 2 // 2 , then, 

(0K|7r b ) T =(l- <0K|7r 6 ) T=0 . (13) 

Since the matrix element is lorentz covariant at zero temperature, it remains 
so to ^T 2 // 2 , and/* = /«. 



"Strictly speaking, the Gibbs average is only defined for diagonal matrix elements but, to 
leading order in T 2 at least, (JlOf) reproduces the known answer. It would be interesting to 
know what is the definition valid to any order. 
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3 Nonlinear sigma model 

It would of great interest to perform calculations in the nonlinear sigma model, 
to ask: how does the damping rate vanish about zero momentum? and, is 
v < 1? Both effects can first appear at two loop order, ~ T 4 j ll ._ .Galculat ions 
of the damping rate have been performed by Gavin et al. ; l!3llJ't3 however, 
they used various (physically reasonable) approximations to the tttt scattering 
amplitude, and so it is not apparent how "f(p) vanishes as p — > in the chiral 
limit [| 

We ,nate, however, that there is indirect support for v < 1 in the results of 
Schenk.tiHl3 He also computed not in the strict chiral limit, but using various 
approximations to the 7T7T scattering amplitude. In fig. 5ofRef. [11] and fig. 7 
of Ref. [12], Schenk graphs f(p) = uj(p)/y /, p 2 + (rn°) 2 , where roj is the mass 
at zero temperature, and u>(p) his computed quasiparticle pion mass. 

At leading order, w; o (0) > m°, and so /; o (0) > 1. As p increases, 
fioip) decreases monotonically to one. This is what happens when 0Ji o {p) — 
V 'p 2 + Wio(0) 2 , or v to = 1. 

At next to leading order, w n ; o (0) < mj. This is not particularly suprising. 
However, as p increases, f n io(p) first decreases, and then increases, approaching 
one from below. That is, there is a "dip" in f n io(p)- Assume that the dispersion 
relation is 

UJnloip) 2 = V n l (p) 2 p 2 +UJnlo(0) 2 ■ (14) 

Of course v n i (p) — > 1 for p ^> w„; o (0). For the dip to occur, however, at zero 
momentum the velocity must satisfy 

Vnlo(0) < fj < 1 • (15) 

The first inequality guarantees that a dip occurs in f n io(p)', the second is 
because Schenk's results show that ui nlo (0) < mj. These are both strict in- 
equalities, and so w„Zo(0) < 1. This is in accord with our analysis. 

4 Linear sigma model 

In this section we describe, hopefully in a more comprehensible way than in Ref. 
[1], where these effects first show up in a linear sigma model in an expansion 
about zero temperature. The modifications o£_pion propagation were already 
computed long ago by Itoyama and Mueller; til the only thing which we are 
doing is to turn their results for the pion propagator into results for /* and 

'To leading order in a virial expansion, as eq. (2.4) of ref. [11], we estimate 7 ~ p(T 4 //*) 
about zero momentum in the chiral limit. 
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The lagrangian of the linear sigma model is given by 



C=\{d,cf>f + \{^f - ha , (16) 

where <f> — (c> ??) is an 0(4) isovector field. We introduce a background mag- 
netic field h which is proportional to the current quark mass m. For h = 0, 
the vacuum expectation value of the a is <jq — y/ /i 2 /A, where we then shift 
cr — > (To + cr; for two flavors, = Co- 

We want to compute terms of order ~ T 4 . There are many such terms, 
both at one and at two loop order. To avoid being drowned in calculational 
details, we extract a small subset of terms which are dominant in weak cou- 
pling, those ~ T 4 /(f 2 m 2 ) at one loop order. In weak coupling, m 2 = 2A/ 2 , 
and so these terms are larger by 1/A than terms ~ T 4 // 4 . Indeed, we are 
lucky not to have to go to two loop order in the first place, as is required in 
the nonlinear model. ._. 

We first review the calculations of Itoyama and Mueller, til There are sev- 
eral diagrams which contribute to the pion self energy, fig. 5 of Ref. [13]. 
However, the only one which can produce the effects we are looking for is 
where a pion splits into a pion and a sigma, and then recombines. This is be- 
cause we need a nontrivial momentum dependence in the pion self energy, and 
this is the only diagram that has it. The integral for this particular diagram 
is proportional to 

where tr K = T^=-oo / d 3 k/(2nf. 

For the real parts of the pion self energy, the trick is simply to expand 
this integral in powers of the external momentum P. At lowest order, this is 
trivial, since we just take the propagator for the sigma meson, and replace it 
by a constant: 

1 1 IT 2 

I P „ tr K — ~ — — . (18) 
m% K m% 12 

We can forget about any ultraviolet divergences along the way, since they will 
be taken care of by zero temperature renormalization, as usual. (We might 
have to worry if we found any logarithmic ultraviolet divergences, but we're 
only throwing away powers divergences, which is ok.) 

To go beyond leading order, then, to ~ T 4 , all we have to do is to expand 
the integral in P to ~ P 2 . Doing so, we find that there are two types of 
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integrals which enter. One is that of Eq. (|l8|), while the other is 

K^K V 7r 2 T 4 
trK ^0~ ~ - ^n") -95- - ( 19 ) 

where n M = (1,0). This appearance of the timclikc vector is entirely where 
v < c will enter at this order, since it produces different terms in the pion self 
energy ~ p\ and p 2 . For both integrals in Eqs. ( |isj ) and ([191) , evidently the 
scale of the loop momentum K ~ T. This justifies the approximations made. 

These are all the integrals which are needed to get the real part of the pion 
self energy. For the imaginary part, we need the integral 

,2 



1 f m'i\ 



Im T(P)U P<m<7 ~ — exp [ ) . (20) 



The imaginary part is exponentially small because one needs to scatter off of a 
sigma meson in the thermal distribution at low temperature. Sigmas are heavy 
at low temperature, so the probability to pop one out of the distribution is 
Boltzmann suppressed. 

Introducing the quantities 

T 2 ir 2 T 4 1 m% ( ml \ 

* x = To72 ' t2 = 7772-^" ' *3 = r^- 72-2 exp I -7-= , 21 
12/ 2 45 f£mi 32tt/ 2 p 2 \ 4 P T / 



to this order the inverse pion propagator is 

Z n A~ x (P) ~ (l + t 1 + 6t 2 )p 2 + {l + t 1 -2t 2 )p 2 +ml (1 + 3U/2) ~2ip% . (22) 

The zero of this expression determines the pion mass shell. In the chiral limit, 
it is 

ip° ~ vp - ipt 3 , (23) 
where the pion velocity is given by 

v 2 ~ 1 - 8t 2 , (24) 

The result for the pion velocity agrees with (A. 8) of Ref. [13]. There appears 
to be a misprint for the pion damping rate in (A. 8') of Ref. [13]. 

To evaluate f* and in the linear sigma model we need the axial current, 

Al = (oo + <7)^7T° - 7T a d^ . (25) 

Then we take the two point function of axial currents and saturate it with a 
single pion in the intermediate state; from Eqs. (9) and (11) of Ref. [1], this 
equals 

(O|0M»0"4|O)t ~ 5 ab {Refl) (fip 2 + Hp 2 ) , (26) 



Generally, for small momenta this two point function is dominated by the 
contribution of a single pion, since the contribution of other states is of higher 
order in momenta, ~ (p 2 ) 2 , etc. This remains true even if the condition of Eq. 
(0) is not satisfied: even for strongly damped pions, the dominant contribution 
to Eq. (|2tj ) arises from single pions in the intermediate state. 

At one loop order, the diagrams which contribute to j\ and /* are given 
in fig. (5) of Ref. [15]. As for the pion self energy, the only diagrams which can 
produce the effects which we are looking for are those involving a virtual pion- 
sigma meson loop. In the calculation of the two point function in Eq. (p6|), 
this kind of loop enters in two ways. First, simply as self energy insertions on 
the intermediate pion. Secondly, from the form of the axial current, they can 
also enter through either vertex, as a type of form factor. 



For the real parts of the pion decay constants, the integrals in Eqs. (18) 
and ( |l9| ) suffice. For the form factor diagram involving a loop with a sigma 
meson and a pion, again one expands the propagator of the sigma meson in 
powers of the external momentum. As a form factor, only terms linear in P 
need be retained. 

For the imaginary parts, we need two new integrals, 

trK K^{P-Kf + m l) ~ lk (if + T ) CXP {-^f) ■ (2?) 



and 

trK K*((P-K)2 + ml) ~W^{^~ T ) CXP [~*pf ) ■ (28) 

These expressions are valid near the light cone, uj ~ p <C m a . In these integrals, 
the dominant loop momentum k is very large, k ~ rri^/p 3> m a . Thus one 
typically retains only those terms ~ k and not ~ p. 
The results for /* and are 

/* ~(l-ti + 3t 2 + its)/* , #~ (1 - h - 5t 2 - it 3 ) U ■ (29) 

It is then simple to verify that Eqs. (||) and (JsJ) give the same mass shell 
condition as found directly from the pion propagator in Eq. (p3[). 
Given that 

(qq) - a (T) - a (0) (1 - 3*i/2) , (30) 

we also verify our generalization of the formula of Gell-Mann, Oakes, and 
Renner in Eq. (||) for the dynamic pion mass: 

m 2 (T) ~m 2 (l + ii/2-6t 2 ) . (31) 
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From the propagator, this result is really trivial, since it follows simply by 
using the wave function renormalization of the pion, Z v ~ 1 + t\ + %ti, and 
finding the pole in p^. In detail, what matters is that the coefficient of ti in 
Z^ and in (fl) 2 are equal. 
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